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ABSTRACT 

The light from shock breakout of stellar explosions carries a wealth of in- 
formation on the progenitor and the explosion. The observed breakout signa- 
ture strongly depends on the velocity of the shock at the time of breakout. 
The emission from Newtonian breakouts, typical in regular core-collapse su- 
pernovae (SNe), was explored extensively. However, a large variety of explo- 
sions result in mildly or ultra relativistic breakouts, where the observed sig- 
nature is unknown. Here we calculate the luminosity and spectrum produced 
by relativistic breakouts. In order to do so we improve analytic description of 
relativistic radiation mediated shocks and follow the system from the break- 
out itself, through the planar phase and into the spherical phase. We limit 
our calculation to cases where the post breakout acceleration of the gas ends 
during the planar phase (i.e., the final gas Lorentz factor < 30). We find 
that spherical relativistic breakouts produce a flash of gamma-rays with en- 
ergy, Eb , temperature, Tb Q , and duration, t^ s , that provide the breakout radius 
(» 5^ (t^/lO s) {T bo /50 keV) 2 ) and Lorentz factor (« T bo /50 keV) and that 
always satisfy (t° b b o s /20 s) ~ (E bo /10 46 erg ) 1/2 (T bo /50 keV) -2 ' 68 . The breakout 
flare is typically followed, on longer time scales, by X-rays that carry a compa- 
rable energy. We apply our model to a variety of explosions, including la and .la 
SNe, AIC, energetic SNe and gamma-ray bursts (GRBs). We find that all these 
events produce detectable gamma-ray signals, some of which may have already 
been seen. Some particular examples are: (i) Relativistic shock breakouts provide 
a natural explanation to the energy, temperature and time scales, as well as many 
other observed features, of all low luminosity GRBs. (ii) Nearby broad-line Ib/c 
(like SN 2002ap) may produce a detectable 7-ray signal, (iii) Galactic la SNe 
produce detectable 7-ray flares, if their progenitors are single degenerate. We 
conclude that relativistic shock breakout is a generic process for the production 
of gamma-ray flares, which opens a new window for the study and detection of 
a variety of stellar explosions. 
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Introduction 



The first observable light from any stellar explosion is generated by the breakout of the 
shock that traverse the star. In most stellar objects the density drops sharply near the stellar 
edge. The shock that encounters this density drop accelerates, leaving most of its energy 
behind. Therefore, it breaks out at a velocity that is significantly higher than the typical 
shock velocity in the stellar interior, carrying only a small fraction of the total explosion 
energy. In typical core-collapse supernovae (SNe), where the total explosion energy is ~ 10 51 
erg and the progenitor is not a compact object, the velocity of the shock at breakout is 
Newtonian. For example in a red supergiant progenitor the breakout velocity is in the range 
5, 000 - 10, 000 km/s, while in a Wolf-Rayet progenitor it is ~ 30, 000 - 100, 000 km/s. 
In cases where the explosion is more energetic or the progenitor is a compact object, the 
breakout can become mildly or even ultra relativistic. 

Before its breakout, the shock is dominated by radiation. During Newtonian breakout 
all the photons that are contained in the shock transition layer are escaping the system 
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properties depend mostly on breakout radius and on the shock velocity at the layer that 
dominates the breakout emission, vq. Among these properties the observed typical photon 
frequency, denoted here as the radiation temperature T , can be especially sensitive to the 
shock velocity. In relatively slow shocks, vq < 15,000 km/s, the radiation in the shock 
breakout layer is in thermal equilibrium and T n oc vl/ 2 : produ cing typically a spectrum that 
peaks in the UV. However, as was shown by IWeaverl (119761 ). in faster radiation mediated 
shocks the radiation is away from thermal equilibrium in the immediate shock downstream, 
and as a result the temperature in the shock transition layer is much higher than the one 
obtained assuming thermal equilibrium. As a result the breakout in fast Newtonian shocks, 
v n ~ 30, 000 - 100, 000 km/s, peaks in X-rays (katz. Budnik k Waxmanlboiol : iNakar k Sari 



20101 ) . 



The emission following the breakout can be divided into two dynamically different 
regimes. At first, before the expanding gas doubles its radius, the dynamics is planar while 



at later time the evolution is spherical (IPiro. Chang k Weinberg |2010| ; INakar k Saril 12010 



Katz. Sapir k Waxman! 12011). The evolution of t he temperature and luminosity during the 



planar phase was studied by lNakar k Saril (120101 ). were we have shown that if the breakout 
emission is out of thermal equilibrium at breakout, it remains out of equilibrium also dur- 
ing the planar phase and thermal equilibrium is achieved only during the spherical phase. 
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All these studies considered only explosions where the shock breakout velocities are Newto- 
nian and pairs are not significantly produced in the transition layer, i.e., (3 = vq/c < 0.5 
where c is the speed of light. 

A large range of cosmic explosions are expected to produce breakout velocities that are 
either mildly relativistic o r ultra relativistic. A few exam ples are a range of white d warf 



explosions such as la SN (IPiro. Chang fc Weinberg 
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2010 



2006 
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The energy involved in most of these explosions is < 10 51 erg, but the compact progenitor and 
relatively low ejecta mass results in relativistic breakouts. Another example is extremely en- 
ergetic explosions of massive stars, with explosion energies in the range 10 52 — 10 53 erg. These 



were recently observed in g rowing numbers, e.g., SN2005ap and SN2007bi (jGal-Yam et al. 



20091 ; iQuimby et al. 1 12009k e.g.,), and are believed to be the outcomes of various progeni- 
tor systems. Yet another example is massive star explosions that are associated with long 
gamma-ray bursts, where the breakout is ultra-relativistic. Relativistic breakouts may also 
take place in choked GRBs, where the relativistic jets fail to penetrate through the stellar 
envelope, and no regular long GRB is produced. As we show here, this mechanism is a 
promising source of all the observed low-luminosity GRBs. 

The physical processes involved in producing the observed signature of a relativistic 
shock breakout is different from the Newtonian case in two fundamental ways. First, once 
the velocity of a radiation mediated shock becomes large enough, /3 B > 0.5, and the tem- 
perature behind the shock exceeds 50 keV, pair production becomes the dominant source of 
leptons, which in turn are the main source of photons via free-free emission. The exponential 
sensitivity of the number of pairs to the temperature regulates the rest frame temperature in 
the shock immediate downst ream to an almost constant value of 100 — 200 keV, in dependent 
of the shock Lorentz factor (IKatz. Budnik &: Waxman! |2010| ; iBudnik et al. Il2010i ). The gas 
becomes loaded with pairs that dominate its optical depth, resulting in a temperature de- 
pendent opacity per unit of mass. This is in sharp contrast to Newtonian shocks, especially 
when 0.05 < f3 s < 0.5, where the temperature of the shocked fluid depends strongly on 
the shock velocity but the opacity per unit of mass is constant. The second major differ- 
ence of relativistic explosions is the hydrodynamics both before and after the breakout. A 
Newtonian shock that pro pagates in a power-law decreasing density assumes the self similar 
solution of ISakurail (Il960l ) before breakout, and the gas velocity is roughly constant after 
the breakout (it doubles by about a factor of two). Once the shock becomes relativistic it 
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follows the soluti on of iJohnson fc McKed (ll97ll ) (see also iTan. Matzner fc McKed l200ll and 
Pan fc Sarill2006l ) before the breakout and each fluid element continues to accelerate signif- 
icantly following the breakout. These two differences, together with the usual relativistic 
effects such as relativistic beaming, govern the observed emission, which as we show here is 
very different than in the Newtonian case. 

Yet another difference between Newtonian and relativistic radiation mediated shocks is 
the physical width of the shock. In radiation mediated shocks t s ~ c/v s , where r s is the total 
optical depth (including pairs and Klein-Nishina effects) seen by a photon going from the 
shock downstream towards the shock upstream. In Newtonian radiation mediated shocks, 
where no pairs are produced, the optical depth of a layer remains similar before it is shocked 
and once it becomes the shock transition layer. Therefore, the shock breaks out at the layer 
in which r = r s , where r is the Thompson optical depth to the stellar edge before the shock 
crossing (i.e., without pairs). However, in relativistic shocks, where pairs are produced 
and Klein-Nishina effects may become important, the optical depth of a layer is changed 
significantly when it is swept up by the shock and the relation between t s and r is not trivial 
anymore. In order to find this relation one needs to know the structure of the shock transition 
layer. The struc ture of relativi s tic rad iatio n mediated shocks in di ff erent regimes was solved 
numerically by IBudnik et al. I (120101 ) and iLevinson fc Bromberd (120081 ). The solution of 
Budnik et alTTfhoioh . where a significant number of pho tons is generated within the shock. 



is the r elevant one to the typ e of shocks we consider here (IBromberg. Mikolitzky &: Levinson 
20111 ). IBudnik et al. I ( 20101 ) also provide an approximate analytic description of the shock 



structure as function of the optical depth. We derive a more accurate analytic description of 
the structure of the shock transition layer and use it to find the value of r at the point that 
the shock breaks out. We find that while production of pairs results in a shock that breaks 
out of the star at t < 1, the observed emission is always dominated by the layer where 
t « 1. We therefore use the subscript V and the terminology 'breakout shell' to denote the 
t = 1 layer, which is not the actual layer where the breakout takes place, but is the layer 
that dominates the observed breakout emission. 

In this paper we calculate the evolution of the observed luminosity and temperature 
from explosions in which the shock becomes mildly or ultra relativistic, i.e., 7oA) > 0.5. We 
follow the light curve as long as the observed radiation is generated by gas that is moving at 
relativistic velocities. At lat er times the radiatio n is determined by Newtonian gas, which 
light curve was discussed in iNakar fc Saril (120101 ) . Our solution is limited to cases where 
the breakout shell ends its post shock acceleration before it doubles its radius. This limit is 
translated to a final (post-acceleration) Lorentz factor of the breakout shell < 30. We also 
assume that the opacity of the progenitor wind is negligible and do not affect the breakout 
emission. Our calculations are applicable to a wide range of energetic and/or compact 
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explosions including la and .la SNe, energetic core collapse and pair instability SNe and 
possibly part of the phase space of chocked GRB jets. 

We find that typically, relativistic breakouts produce flashes of gamma-rays that with- 
out further information (e.g., the burst redshi ft) can be mist aken as cosmological gamma- 
ray bursts. These results provide a closure to IColgatd (119681 ). who was the first to predict 
that shock breakouts produc e flares of gamma-rays, even befo re the detec tion of g a mma- 
ray bursters was announced (IKlebesadel. Strong fc Olsonl 119731 ). However, IColgatd (119681 ) 
assumed that the gas behind the shock is in thermal equilibrium, and the observed 7-rays 
are due to blue shift of ultra-relativistic gas, brought to high Lorentz factors by shock 
that accelerates indefinitely as it approaches the stellar surface. It ignores deviation from 
thermal equilibrium and the fact that acceleration stops once the shock width is compa- 
rable to the distance to the stellar surface. As it turned out, the deviation from ther- 
mal equilibrium increases the photon's rest frame temperature and compensates for the 
limited Lorentz factor of the gas, resulting in a burst of 7-rays. Since Colgate's predic- 
tion, breakouts were suggested to pro duce bright gamma-ray emission by several authors 



in the context of low-luminosity GRBs ( Kulk arni et al. 



Campana et al. 1120061 ; IWaxman. Meszaros fc Campana 
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Wang et al. 1120071 ). All these 



calculations assume thermal equ ilibrium, underestimating the true radiation temperature. 
Katz. Budnik fc Waxmanl (120101 ) realized that the deviation from equilibrium, discussed by 
Weaver! (119761 ). will lead to 7-ray breakouts. They did not follow, however, the post shock 
dynamics and opacity evolution to find the properties of the observed emission. Here we 
provide, for the first time, a quantitative description of the predicted light curve from rela- 
tivistic breakouts. We show that these light curves explain very well many of the observed 
features of low-luminosity GRBs, and also, that a variety of other stellar explosions are most 
likely producing a bright gamma-ray breakout emission. 

In §2] we discuss the pre and post breakout dynamics. The dynamics is dominated 
by mass, energy and momentum conservation and is therefore largely independent of the 
photon number and spectrum, as long as the gas is optically thick. In £J3] we calculate the 
temperature and opacity evolution of the expanding gas. §4] follows the observed light curve 
from the breakout to the spherical phase. Our results are applied in §S]to find the observed 
signature of known and hypothesized explosions. An analytic description of the structure of 
the transition layer in relativistic radiation mediated shocks is presented in Appendix A. 
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2. Dynamics 

2.1. Pre-breakout dynamics 

Consider a relativistic blast wave that propagates near a stellar edge pre-explosion mass 
density gradient of the form p oc z n , where z = (R* — r), R* is the stellar radius and r 
is the distance from the star center. In compact stars, where shocks are more likely to 
become relativistic, the envelope is radiative and typically n 3, which is the value that 
we use throughout the paper. The shock acceleration stops when the width of the shock 
is comparable to z and the shock breaks out. As we show in appendix A, this takes place 
in the shell where the pre-explosion optical depth to the stellar edge is r ~ 0.0l7 s , where 
7 S is the shock Lorentz factor. Namely, breakout takes place at z ~ 0.01^ s /(pkt), where 
Krp ps 0.2 cm 2 /g is the Thompson cross section per unit of mass^- This sets the maximal 
Lorentz factor of the shock and we refer to that shell as the outermost shell. Unlike the 
Newtonian case, the observed emission from relativistic breakout is not dominated by the 
outermost shell. Instead it is dominated by the shell where the pre-explosion optical depth 
ia^| r ~ 1 (i.e., z ~ l/(pfi>r)), which we therefore refer to as the breakout shell and denote 
its properties by the subscript V. In case of properties that evolve with time, this subscript 
refers to their value in the breakout shell right after it is crossed by the shock. Given 
the above assumptions the system is completely defined by the following parameters of the 
breakout shell: 

• 70: the Lorentz factor of the shock when it crosses the breakout shell. 

• m§. the mass of the breakout shell. 

• d Q : the post-shock width (i.e., after compression) of the breakout shell in the lab frame. 

• R*\ the radius of the star. 

• n: the power law index describing the pre-explosion density profile, here we use n — 3. 



1 Throughout the paper we use this value, which is appropriate for hydrogen free ionized gas. The 
dependence on the value of kt is weak and all our results are applicable also to hydrogen rich ionized gas. 
Note that through most of the phases discussed here the temperatures in the system are larger than 10 
keV, where the gas is fully ionized regardless of its composition, implying that the results are insensitive to 
mctallicity. 

2 This is the criterion for the dominating shell as long as 70 J$ 100. As we discuss below this criterion is 
always satisfied in the regime at which our analysis is applicable 
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Other characteristics of the breakout shell can be calculated from the above. For example, 
the initial thermal energy of the breakout shell is E ps 7 2 moC 2 . 

The pre-breakout and post-breakout dynamics of a planar relati v istic blast-wave in a de- 
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similar both before and after the breakout. Before breakout a relativistic shock accelerates 
as 7j oc where p = (\/3 — |) ps 0.23 and the subscript V indicates initial values (note 
that p is the pre-shock density). This propagation profile sets the entire pre-explosion hy- 
drodynamic profile. Conceptually, it is useful to treat the whole expanding envelope as a 
series of successive shells. The properties of any shell deeper than the outermost shell, can 
be characterized by its mass, m: 

• 7» = 70(^/^0) "+ 1 oc m -0 ' 17 is the initial Lorentz factor. 

• di ~ z/7 2 oc m (" +1 > ps m - 6 is the width of a shocked shell at shock breakout, in the 
lab frame. 

l-\-nfl n Art 

• d\ — d^i oc m< n+1 ) ~ m is the width of a shocked shell at shock breakout, in the 
shell rest frame. 

• Ei ps 7?mc 2 oc m "+ 1 ps m - 65 is the internal energy of a shell at shock breakout, in 
the lab frame. 

• ti ps ^-7^ oc m"+i oc m ' 25 is the lab frame time for the shell expansion, which is also 
the time between the shock crossing and the breakout 



2.2. Post-breakout dynamics 

Following breakout all the shells expand and accelerate. The acceleration is faster than 
that of a freely expanding relativistic shell as outer shells gain energy from inner ones. In 
this subsection we follow the acceleration of shells with mass m > mo- These shells remain 
opaque during the whole planar phase, also after their pairs annihilate. The final Lorentz 
factor of less massive shells depends on their evolving opacity and is discussed shortly in the 
next section. 



The acceleration during the planar ph ase, i.e., before a shell doubles its radius, follows 
(Ijohnson fc McKee!E97ll IPan fc Sarilhood ): 
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where t is the lab frame time (not to be confused with observer time) measured since break- 
out. If the shell is optically thick at the end of acceleration and acceleration ends during the 
planar phase then the final Lorentz factor is 

7/ = oc m- ' 48 . (2) 

This relation is general as it does not depend on the exact density profile, as long as there 
is a large energy reservoir behind the accelerating shell. The acceleration ends at 

t f = taf +V * oc m" ' 57 . (3) 

Thus, more massive shells end their acceleration at earlier times and lower Lorentz factors. 

In this paper we restrict our treatment to cases where the breakout shell ends its accel- 
eration during the planar phase, i.e., tfp < t s where 

*. = - (4) 
c 

is the lab frame time of transition between the planar and spherical phases. Since pre-shocked 
optical depth of the breakout shell is of order unity (i.e., KtPqZq ~ 1) : 



h = £1 ~ ( J*L) n+1 = 3 . lO^M^ - 25 ^- 5 
t a R* \k t MJ 5 5 



(5) 



where M x and R x are the progenitor star mass and radius in units of x solar masses and 
radii respectively. Thus, there is a critical Lorentz factor, 7o jS , for which tfp = t s \ 

7o ,s = 3.5 M^Rf-\ (6) 

The corresponding critical final Lorentz factor of the breakout shell is: 

7/,. = lfs +l = 30 M^Rf * (7) 

If 7o < 7o,s then acceleration ends in the planar phase and the corresponding final Lorentz 
factor is given by equation |5J Otherwise it ends during the spherical phase and equation [2] 
is not valid anymore, a case that we treat elsewhere. 

Acceleration continues beyond the breakout shell, up to the outermost shell, which 
satisfies z ~ 0.01^/ (p^t)- This shell achieves the maximal initial Lorentz factor: 

lma X ,i « 2 7o °- 85 (8) 

Note that our calculation is limited to 7 < 7o, s <C 100, where 7 < 7 ma z,i and the entire 
above description is consistent. 
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3. Temperature and opacity 

The r est frame temperature behind a radiation mediated sh ock with /3 S > 0.5 is roughly 
constant fcatz. Budnik fc Waxmanlhoiol : budnik et oTIboioh : 



T[ ~ 200keV (9) 

The reason is that at these velocities the radiation behind the shock is out of thermal 
equilibrium and is set by t he ability to g enerate photons. The main photon source in the 



shock is free- free emission (I Weaver! 1 19761 ). Once the temperature behind the shock rises 
above ~ 50 keV, pairs are generated and their number becomes dominant over the number 
of electrons that are advected from the upstream. These pairs significantly increase the rate 
of photon generation. The strong dependence of the number of pairs on the temperature 
in that range serves as a "thermostat" that sets the temperature in the immediate shock 
downstream to be nearly constant, regardless of the shock Lorentz factor. Although the 
radiation is out of thermal equilibrium the photons and the pair loaded gas in the immediate 
downstream are in C ompton-pair equilibrium and there is a single gas-radiation temperature 



( IBudnik et al. II2010I ). This property of relativistic radiation mediated shocks is dominating 
the emission at the shock breakout, during the planar phase and the relativistic part of the 
spherical phase. 

The gas behind the shock is pair loaded and therefore its optical depth is very high, 
preventing any significant leakage of photons, even if the pre-shocked (unloaded) gas optical 
depth is below unity (the pairs set a diffusion time that is much larger than the dynamical 
time). On the other hand, the rest frame temperature in the expanding gas drops, and with 
it the number of pairs and the photon generation rate. Therefore, the number of photons in 
the opaque expanding shells is roughly constant. The mildly relativistic temperature and the 
large optical depth ensure that photons and pairs interact many times (i.e., sharing energy, 
annihilation and creation) over the dynamical time, thereby keeping the pairs-photon gas in 
Compton pair equilibrium. 

The radiation of the breakout shell is confined to the gas during its expansion until 
the rest frame temperature drops enough so the pair loading becomes negligible. During 
the expansion of a shell its rest frame temperature falls as T' oc V' 13 , where V is the shell 
volume (in its rest frame) and 1 + (3 is an effective adiabatic index. The value of ft depends on 
temperature and it can drop slightly (by up to 30%) below 1/3 because of pairs production 
and annihilation. Here we neglect this small deviation from the typical relativistic equation 



3 This is true as long as the blackbody temperature that corresponds to the post shock energy density is 
lower than T[ 
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of state and approximate /3 = 1/3. During the planar phase the volume of a fluid element 
grows as V oct/y. Therefore during acceleration V oc t( 3_v/3 ")/ 2 while after the acceleration 
ends V oc t. During the spherical phase V oc t 3 implying: 



T' = T'l 
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t, < t 



(10) 



The optical depth of the breakout shell drops to unity once the pairs density, n±, is 
lower than the proton density, n p (and their accompanied electrons ), i.e., n±/n p < 1. The 
initial value of this ratio is (e.g., ISvenssonlll984[ iBudnik et al. 1120100 : 

„2 



n±,o 7o™ p <r 



10 J 7o 



11) 



n pfi 3kT( 

Since pairs are in annihilation-creation equilibrium their density drops exponentially with 
T' at this temperature range and the shell opacity is dominated by the electrons that were 
advected from the upstream once its temperature is: 



1 th ~ 



50 keV, 



(12) 



where the dependence on 70 is very weak. Around this temperature, the breakout shell, and 
any less massive shell, becomes transparent and its radiation escapes. We define t t h,o as the 
time that the breakout shell becomes optically thin, i.e., its T' = T[ h . 

There are three critical times in the system, in which the dynamics and radiation change 
their behavior. These are the acceleration and transparency times of the breakout shell, tf$ 
and tth,o respectively, and the planar to spherical transition time, t s . The observed signal 
depends on the relative order of these three time scales. Here we restrict our analysis to the 
case that acceleration ends during the planar phase, i.e. 
equation [10] implies: 

S + V^ rpl 



tffi < t s . Under this assumption, 
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(tf,o < tth,a < ts) 
(t s < tth,o) 



(13) 



Now, Tq/T^ ~ 4 and we consider only 70 < 70, s , which for any non-degenerate progenitor 
implies 70 < 7o, s < 4. Therefore, the breakout shell becomes transparent only after accelera- 
tion ends, i.e., tf )0 < t t u,o- On the other hand, the opacity drops quickly during the spherical 
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phase, so equation [131 implies that the breakout shell becomes transparent during the planar 
phase, or at the beginning of the spherical phase, namely t t h,o i$ ts- Therefore for a large 
range of realistic explosion parameters tf t0 < t th < t s , which is the case that we present in 
detail below. Finally, note that shells with 7^ > 70 can be accelerated only as long as they are 
optically thick, namely their T' > T[ h . Therefore, only shells with 7, < min{4, 7o,s} satisfy 
the relation 7/ = 7^^, shells with higher 7, end their acceleration before this relation is 
satisfied. 



4. Light curve 

4.1. Breakout emission 

The radiation in the expanding gas is trapped by pairs at the time of breakout. In the 
previous section we have shown that for most progenitors in the regime that we consider, 
the breakout shell becomes transparent after it ended acceleration and before, or soon after, 
the transition to the spherical phase, i.e., £/ ; o < Uh,o ^ ts- Therefore we derive below the 
observed light curve in that case. Note that the breakout shell is not the only shell that 
becomes transparent during the planar phase. All the lighter shells (m < mo) out to the 
the outermost shell, become transparent during the planar phase as well . Shells that are 
deeper (and more massive and slower) than the breakout shell remain opaque during the 
planar phase and become transparent only during the spherical phase. 

Since slower moving shells carry significantly more energy, Ei oc 7~ 3 ' 75 , the breakout 
emission is dominated by the energy confined to the breakout shell at tfh,o, when its local 
frame temperature is T' th : 

E bo = E ^^ » \m c 2 l0lf « 2 • l^Rhp 1 , (14) 
7o J-o 4 

where we use r ~ KT^o/(4vri?^) = 1 to find 

m w 4 • lQ- 9 M Q Rl . (15) 
The observed temperature of this emission is 

T bo = T^ffi ~ 507/ keV , (16) 
and it is smeared by light travel time effects over a duration that is comparable to t° bs : 

C « ^^h^lO^s, (17) 

c 7/,o 7/, 
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where t° bs is the time of transition from the planar to the spherical phase as seen by the 
observer. Clearly, equations [ToTTlTl are general and are applicable to any gas density profile. 
In fact also equation [TJ] is applicable to any sharply decreasing density gradient, even if it 
is not a power-law. The reason is that m does not depend on this profile and so does the 
relation between 70 and 7^0 (equation [2]), as long as acceleration ends during the planar 
phase and there is a large energy reservoir behind the breakout shell. Together, equations 
[1411171 provide three generic observables that depend on only two parameters 7^0 and R*. 
Thus, the energy, temperature and time scale of the breakout flare are enough to determine 
-R* and 7^0, and they also must satisfy 

1/2 9+V3 

* ~ 20s (lo^) («?Sv)~ * • < 18 > 

if the source of the flare is a quasi-sphericafl relativistic shock breakout. This relativistic 
breakout relation can be used to test the origin of observed gamma-ray flares. 

During t < t^ s the observer receives also photons from shells that are outer of the 
breakout shell, which have higher observed temperature than T bo . These photons generate 
a high energy power- law at frequencies T bo < v. Shells with 7^ < T[jT[ h « 4 obtain their 
terminal Lorentz factor before they become optically thin, implying that 7/ = / y/ 3+1 and 
E oc Eijf/ji oc 7~ 2 ' 01 while T = T^ h jf oc / y/^ +1 - Therefore the integrated spectrum of the 
breakout flare exhibit a high energy power-law 

vF v (v > T bo ) oc z/-°- 74 (19) 

up to v = min{2, 0.257q 05 } MeV, which can be an order of magnitude larger than T bo . 
During t Q b s < ^o S (= t° bs ) the observer receives also photons from the break out shell itself 



after it becomes transparent. During this phase adiabatic cooling dictates (e.g-. lNakar fc Sari 



2010b, T oc tJJ 3 and L oc tj/ 3 , implying 

vF v {v < T bo ) oc v (20) 

The range of the low energy power-law in rather limited T bo > v > T bo {ct t hfl/ R*) 1 ^ 3 = 
Q.bTbo^y^ M^ 1 ! 12 R l J Q . Note that equation [19] and the spectral range where the high and 
low energy power-laws are observed, do depend on the specific density profile and are given 
for n = 3. 



4 Ebo, Tbo, and fj^f depend only on emission from a region within an angle of ~ l/7/,o with respect to the 
line of sight 
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At t° bs the observed luminosity is still dominated by light emitted when the breakout 
shell become transparent, ~ ttho- This emission fades quickly until the emission from the 
spherical phase becomes dominant. The fading light curve can be easily derive in case that 
7/ o ^ 1 ; since then the observed luminosity is dominated by emission of the breakout shell 
at large ang les (> l/7/ t o). The luminosity decays then as L oc t~ 2 s and the temperature 
as T oc t (IKumar fc Panaitescul |2000| ; iNakar fc Piranl 120031 ) , until the emission from the 
spherical phase becomes dominant. 



4.2. Spherical phase 

During the spherical phase the radius of the expanding sphere is not constant anymore. 
Instead, the radius of the shells is r oc v t and the luminosity is dominated by photons that 
are leaking from the shell that satisfies r ~ c/v. We re fer to this shell as th e luminosity shell 



sec 



and denote its properties with the superscript ' 

cussion), such that, e.g., the luminosity shell satisfies, by definition, r 



Nakar fc Sari 2010 for a detailed dis 



c/v. The evolution 

during the spherical phase depends on whether the luminosity shell is relativistic or not and 
on the initial temperature of the luminosity shell. As long as 7/ ^> 1 the initial temperature 
of the luminosity shell is ~ 200 keV and its dynamics is well approximated by the relativistic 
self-similar solution. Shells with 0.5 < 7^ < 1 also have an initial temperature of ~ 200 
keV but their dynamics is better described by Ne wtonian approx imation. Namely, shock 
propagation according to the self similar solution of ISakurail (119601 ) and no significant accel- 
eration after shock breakout. The transition time between the relativistic and Newtonian 
phases is calculated below and it takes place around: 

;5* 1.05 _ f obs 3.05 

„ 7/,o - l s 7/,o 



t 



obs 
NW 



(21) 



During the Newtonian phase t f, s oc v~ A . The end of the phase in which the initial tempera- 
ture of the luminosity shell is constant takes place when v ~ 0.5c, which is at t b s w lOffiw- 
Here we present luminosity and temperature evolution during the spherical phase up to that 
point (v ~ 0.5c). At later times the temperature drops quickly until thermal equilibrium 
is obtained. Therefore there is a sharp break in the temperat u re evo lution at 10tAryi/. The 
light curve evolution at t > lOtNw is covered by lNakar fc Saril ( j2010l ). 



4.2.1. 



Relativistic phase (t° bs < t obs < t^yy ) 



The relativistic velocity of the shells dictates r ~ 1 and since r oc m/r 2 the mass of 
the luminosity shell is m oc t 2 . The photons from the luminosity shell arrive to the observer 



-14- 



at t bs ~ t/ij and their arrival is distributed over a comparable duration. Therefore the 
luminosity shell satisfies m oc t°^ 9 , % oc t~ b ° s 12 , jf oc t~ b ° s 33 and £^ oc nvji oc t°^ 7 , while 
f oc t oc t 6s7/ oc t°^ 4 and equation [10] dictates T' oc i^' 3 ■ The observed luminosity and 
temperature are: 



E _ ZffEKT'/T-) _ E c ft y* 14+2 



V3 ^ a / + . \ -1.12 



7o"^ ( ^ ) (22) 



tobs ^obs R* \ts J V 



to A 1/3 * /W- - 68 



T « 200 keV ( ^ I (23) 



The relativistic phase ends once 7, ~ 1 at a time given in equation 



4-2.2. Newtonian phase with constant initial temperature (t°^ < t b s < 10 t^iW 

The luminosity during this phase is independent of the temperature since pairs are long 
gone. It is therefore similar to the luminosity evolution of the spherical phase following a 
Newtonian breakout, where: 

L oc t- b ° s 35 (24) 

Therefore a significant flattening of the luminosity evolution is expected at t^w- From this 
point the luminosity decay is steady until recombination and/or radioactivity start playing 
a role. 

The initial volume of the luminosity shell is proportional to its initial width (all shells 
starts at i?*), rfj oc t 0Ai . The observed temperature evolves as: 

This slope is not very different than that of the relativistic phase. Therefore there is no sig- 
nificant feature in the temperature evolution at t^w Only at t t, s ~ 10 t^w the temperature 
starts dropping rapidly until the luminosity shell gains thermal equilibrium. 



5. Predicted signal of various explosions 



We use the model developed in the previous section to calculate the signal expected 
from known and hypothesized explosions. There is a large number of explosions, such as 
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various SNe types, in which the bulk of the mass is moving in Newtonian velocities but the 
shock accelerates a minute amount of mass to a relativistic breakout. Such explosions are 
often observed mostly in the optical, where the kinetic energy, E, that is carried by the bulk 
of the ejected mass, M ej -, can be measured. Therefore we first provide the characteristic 
properties of the breakout and following emission of such explosions as a function of E, M e j 
and R*. Then we use these results to discuss the predicted signals from various explosions. 

Given E, M e j and the initial Lorentz factor of the breakout shell is determined by 
following the shock acceleration from the Newtonian phase ('jsPs oc p~ 019 ), to the relativistic 
phase (~i s f3 s oc p~ a23 ): 

7o /3o - 2.6 CM-J 45 ^ - 35 (26) 
where E^ = E/(10 52 erg). Th e numerical coefficient is taken to fit the numerical results of 



Tan. Matzner fc McKed (120011 ). which simulate the transition of the shock from Newtonian 
to relativistic. Equation [26] provides a good approximation to the breakout Lorentz factor 
for j3o > 0.5. The initial energy of the breakout shell is 



E « 5 ■ 10 46 erg E^M^R^, (27) 



and its final Lorentz factor is 

7 /,o « UEHM^Rs™ 5 (28) 

where we use the approximation^ 70 ~ 7 oAj- Figure [1] shows a color map of the value of 
7 j o as a function of ejected mass and stellar radius for different values of explosion energies, 
where only the range where our results are valid, 0.5 < 7o A) < 7o,s> is colored. The figure 
also includes the typical phase space location of various explosions in that range. It is 
evident that breakouts of white dwarf explosions, extremely energetic SNe and possibly low 
luminosity GRBs are all within the relevant velocity range. 



5 This approximation is good for 70A) > 2. For lower values of 70 A) it is better to calculate 7^0 = 7o +% ^ = 
(loPo + l) - ^ and plug it into equation [29] 
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Fig. 1. — A color map of the value of 7^0 as function of ejected mass and stellar radius for 
different values of explosion energies. Only cases with 0.5 < 7 /3o < 7o,s are colored. The 
approximate phase space location of various explosions where the breakout velocity falls in 
this range are marked as well. 



Plugging equation EH] into equations [CTTTl we find the observed breakout emission: 



3-^3 

E bo « « 6 • 10 46 ergE^M-^Rr 



T bo « 50 keV 7 /,o ~ 700 keVE^M^R^ 95 



j.obs ^ j.obs * 

Ho ~ H ~ „_,2 



C V,0 



0.06 s E^ A M 2 e -%Rl- 9 



(29) 



L bo « |£ - 4 ■ 10 47 erg/s E^M^^ 
Ho 
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Note that this equation is applicable only to a power-law density profile (n=3) (unlike equa- 
tions [HHT7] which are general), due to the dependence of 7/ j0 on the density profile in case 
that E, M e j and are the physical parameters. 

The breakout emission dominates during the entire planar phase. At the end of the 
planar phase (t° bs ) the luminosity and temperature drops quickly (L oc t~ b 2 s , T oc t~ hs if 
7/ o 3> 1; see discussion in previous section) to join the spherical emission, which is dominated 
by relativistic ejecta until 

tf w « 200 s ElfM-f, (30) 

where we neglect a very weak dependence on R*. The luminosity and temperature at that 
time are 

L{tt w ) « 3 • 10 44 erg^M" / 5 ^- 65 

(31) 

T(tf w ) « 0.2 keVE^M^R™ 

where we ignore very weak dependence^] on M* and use the approximation 7 « 7o/?o- At 
tf s < t obs < t°£ s w the luminosity shell is relativistic and L oc t^ 12 while T oc £~ 6 °" 68 . At 
t°N S w < ^bs < IO^vk the luminosity shell is Newtonian and L oc t~ fe ° s 35 while T oc t~° s e . The 
energy emitted per logarithmic scale in time, Lt Q b s , is rising during that last phase (oc £^ 5 ), 
and most of the energy that is emitted in X-ray (for typical parameters) is emitted around 
t bs ~ lOt^yy- This energy may be comparable to, or even larger than, the breakout energy: 

Lt ob8 {mt w ) « 3 • 10 47 erg E^M^R^ (32) 

At later time the luminosity evolution remains unchanged (until recombination or ra- 
dioactive de cay become domina nt), while the temperature drops rapidly towards thermal 



equilibrium ( INakar fe Saril |2010| ) . 



5.1. White dwarf explosions - la SN, .la SN and AIC 



There are several mechanisms that can result in a shock breakout from a white dwarf. 
The most famous is ty pe la SNe, where t he whole star explodes. Other theoretically predicted 
scenarios are .la SNe (jShen et al. 1120101 ) . where helium accreted in an AM CVn is detonated, 



6 In some explosions M e j — M* while in others M e j < M * . In any case the observations in all the regimes 
we discuss here are almost independent of M* . 
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and accretion induced collapse (AIC) to a n eutron star (e.g., iHillebrandt. Nomoto fc Wolff 
1984 ; iFryer et al. Ill999l ; iDessart et al. 1120061 ) . All these explosions are expected to produce 
a rather constant ejecta velocity {2E / M^) 1 ^ 2 ~ 10,000 km/s, with explosion energy that 



ranges between 10 51 erg in type la SNe, where the whole star is ejected (M e j w 1.4M ), to 
10 49 — 10 50 erg in the other explosions, where only 1 — 10% of the stellar mass is ejected. 
There are also various observed SNe, which are probably generated by a range of white 
dwarf explo sions. Some possib le examples are SN 2005E, where M e j m O.3M and E rj 
4 • 10 50 erg JPerets et a/j kuoj), and SN2010X where M ej « O.16M and E « 1.7 • 10 50 erg 



kasliwal et aTlboioh 



Since E/M e j is similar for all these explosions, and the dependence of 7o/?o o n E (or 
M e j) when this ratio is constant, is very weak (oc E 0A7 ), shock breakouts from all these 
types of white dwarf explosions produce a similar signature. The breakout Lorentz factor 
is 7oA) ~ 1 — 3 for a white dwarf radius 3 — 10 x 10 8 cm. The total energy released in the 
breakout is 10 40 — 10 42 erg and the typical photon energy is ~ MeV. The breakout pulse 
is spread over ~ 1 — 30 ms having a peak luminosity of ~ 10 44 erg/s. These results for 
the breakout emission ar e differ ent, especially in the temperature prediction, than those of 
Piro. Chang fc Weinberg! ( 120101 ) that found a type la breakout signal that peaks in the X- 
rays. The reason is that they assum e thermal equilibriu m behind the shock, which is not 
a good assumption for type la SNe (INakar fc Saril 120101 ). and that they ignore relativistic 
effects. 

Note that this signature is expected only if the line of sight to the observer is transparent. 
In case that the explosion is triggered by white dwarf mergers, the optical depth of the debris 
that surrounds the exploding core can be high and the shock breakout will take place at much 
larger radius than the white dwarf r adius , carrying out much more energy at a much lower 
velocity. For example, IFryer et al. I ( 120101 ) find that a type la SN shock breakout in case of 
a double degenerate merger take place at a radius of ~ 10 13 — 10 14 cm where the shock is 
Newtonian. Thus, a breakout emission from type la SN can potentially differentiate between 
single and double degenerate progenitor systems. 

A flare of 10 41 erg in MeV 7-rays can be detected easily if it takes place in the Milky 
way and possibly also in the Magellanic clouds. Such detection will provide a constrain on 
the on the exact explosion timing as well as on the radius of the exploding white dwarf 
and explosion energy. It may also shed light on the explosion mechanism. This detection 
channel is extinction free and may be the main detection channel in case that the burst take 
place in an obscured Galactic environment. Given the estimated Milky way rate of type la 
SNe (~ 0.01 yr _1 ), the chance to detect such a burst in the coming decade is low, but not 
negligible. The rate of events such as SN 2010X is very hard to constrain observationally, 
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based on current optical surveys. In addition, there are no observational constraints on the 
rate of events with energy lower than 10 50 erg, since their optical signature, where these 
events are searched for, is too weak and evolve too fast for detection. Similarly there are no 
tight theoretical constraints on the rates of .la SNe and AICs. Since the gamma-ray shock 
breakout signatures of all these events are similar and detectable for any galactic event, it is 
worth looking f or them in th e curre nt data of satellites like BATSE, Konus-Wind, Swift and 
Fermi. In fact, ICline et al. I (120051 ) find that the sub-sample of several dozen BATSE and 
Konus-wind GRBs, those with duration shorter than 0.1 s, shows a significant anisotropy 
with concentration of events in the galactic plane (but not towards the Galactic center). This 
sample of events may contain some already observed shock breakouts from white dwarfs. 



5.2. Broad line Ibc SNe 

Broad-line Ibc SNe are thought to be explosions of Wolf-Rayet progenitors, and in some 
cases are very energetic (E ^> 10 51 erg). Some Broad-line Ibc SNe are associated with GRBs 
(these are discussed later), but also those that are not associated with GRBs may produce 
a detectable 7-ray signal. For examp le SN 2002ap has a tota l explosion energy of ~ 4 • 10 51 
erg and an ejected mass of ~ 2.5M (IMaurer &: Mazzalill2010f ). implying a mildly relativistic 



breakout with an energy output of ~ 5 x 10 A4 R® J erg within ~ IOR5 s. At a distance of 
about 7 Mpc the observed ffuence on earth is ~ 10~ 7 i?5 7 erg/cm 2 . This fluence is detectable 
by sensitive gamma-ray detectors even for R* = R Q (although no large field-of-view detector 
with high sensitivity was operational in 2002). Note however, that if a dense wind surrounds 
the progenitor then the flare properties may be very different. The total emitted energy is 
higher while the typical photon energy may fall either within or below the observing band 
of soft 7-ray detectors. 



5.3. Extremely Energetic SNe 

Several types of extremely energetic SNe were detected in recent years. SN 2007bi has 
shown a typical veloc ity of 12,000 km/s and > 5OM of ejecta implying E > 10 53 erg 
(IGal-Yam et al. II2009I ). Most of the optical luminosity observed in this SN is generated by 
radioactive decay. The radius of the progenitor is unknown, but since there is no trace of 
either hydrogen or helium its radius is most likely in the range, R* ~ 10 11 — 10 12 cm. These 
values implies 7o/3o ~ 1- Therefore the breakout is expected to produce a ~ 100 keV flash 
that lasts seconds to tens of seconds and carry an energy of 10 44 — 10 46 erg. 
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Quimby et al. I (120091 ) reports on another type of extremely energetic SNe. They find a 
typical velocity of 14,000 km/s and > 5M of ejecta implying E > 10 52 erg. The observed 
optical energy in these explosions is not generated by radioactive decay and the possible 
sources are interaction of the ejecta with matter that moves more slowly or a central engine. 
In the former case, all the energy is given to the ejecta during the explosion and a bright shock 
breakout is expected. These explosions also do not show any trace of hydrogen and Helium 
implying that their progenitor is also rather compact. Assuming again i?* ~ 10 11 — 10 12 cm 
the signature of the shock breakout from these events is similar to that expected from SN 
2007bi. 

An energy output of 10 44 -10 46 in soft gamma-rays within a minute is detectable out to 
a distance of 3-30 Mpc. Therefore the detection of breakouts from such events is likely only 
if their rate is larger than 10 4 Gpc -3 yr _1 . The rate of extremely energetic SNe is unknown, 
but it is most likely much lower than that since otherwise they should have been detected in 
optical searches in large bunbers. We therefore do not expect current gamma-ray detectors 
to detect any of these SN types, although Nature may surprise us. 



5.4. Low luminosity GRBs 



The engine of long GRBs is thought to be produced during the col lapse of a massive stel- 
l ar cor e and thus to be launching relativistic jets into a stellar envelope. iBromberg. Nakar fc Piran 
(120111 ) show that (at least in some) low-luminosity GRBs it is highly unlikely that the jets 
successfully punch through the star and emerge in order to produce the observed 7-ray 
emission. In that case what can be the source of the observed gamma-rays? 

There are several common features observed in all low luminosity GRBs. All show 
10 48 — 10 50 erg of gamma-rays or hard X-rays that are emitted in sing l e pul sed light curves 



showing a spectra that become softer with time (see iKaneko et al 



2007 and references 



therein). The y also show radio afte r glows that indicate on simila r energy in mildly rela- 



tivistic ejecta jKulkarni et al. I Il998l : ISoderberg et al. I I2OO4L bood ). This energy is only a 



small fraction of the total energy observed in the associated SNe (10 52 — 10 53 erg). The 
smooth light curves and the energy coupled to a mildly relativistic ej ecta motivated severa. 



authors to suggest that the origin of this emission is shock breakouts ( 


Kulkarni et al. 


1998; 


Tan. Matzner & McKee 2001: 


Campana et al. 


20061: 


Waxman, Meszaros & Campana 


2007; 


Wang et al. 


2007J). Later, 


Katz. Budnik & Waxman 


(201oh have shown that the deviation 



from thermal equilibrium can explain the observed 7-rays, compared to the expected X-rays 
when equilibrium is wrongly assumed. However, the question whether shock breakouts are 
indeed the source of low-luminosity GRBs remained open, mostly since there was no quan- 
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titative model that could answer whether shock breakouts can indeed explain even the most 
basic observables, such as energy, temperature and time scales, of these bursts. Below, we 
use our model to answer this question. In section H] we have shown that this three observed 
scales, E bo , T bo and t^ 3 , over-constrain the physical parameters since they depend in the 
case of a quasi-spherical relativistic breakout only on two parameters, 7/^ and the breakout 
radius. Thus, in addition to finding the value of these two physical parameters, the three 
observables must satisfy the relativistic breakout relation (equation [T8|) . This can be used 
as a test that any flare that originates from quasi-spherical relativistic breakout must pass. 
This result is generic as it is independent of the pre-breakout density profile (as long as it is 
steeply decreasing). Below we apply this test to the four observed low-luminosity GRBs 

The four well studied low luminosity GRBs are divided into two pairs with very different 
properties. GRBs 060218 and 100316D are soft and very long. Their ra ther soft spectrum is 
best fitted at early ti me with a cut-off around a peak energy of ~ 40 keV ( iKaneko et al. 1120071 : 
Starling et al. 11201 ll ) and both emitted a total energy of ~ 5 x 10 49 erg. Plugging these values 
into equation [18] we find that if these bursts are flares from quasi-spherical breakouts then 
they should have durations of ~ 1500 s, which is indeed comparable to the observed durations 
of these two bursts. The breakout radius in this model is ~ 5 x 10 13 cm and 7^0 ~ 1 implying 
7oA) ~ 0.3 — 1. In the case of GRB 060218, radio emission suggest breakout velocities on 
the upper side of this range (jSoderberg et al. II2006I ). The large radius inferred by the long 
duration implies that the breakout is most likely not from the edge of the progenitor star, 
but from an opaq ue material thrown by the progenitor to large radii (possibly a wind) prior 



to the explosion (jCampana et al. II2006I ). The fact that a spherical model provides a goo d 



explanation implies that the suggested a-sphericity (IWaxman. Meszaros fc Campanal 120071 ). 
if exist, does not play a major role in determining the burst duration. Note though, that it 
may play an important role in shaping the breakout spectrum above and below T bo and the 
entire light curve following the breakout flare. 

The other two low- luminosity GRBs, 980425 and 031203, are relatively hard, T > 150 
keV, and not particularly long, « 30 s. Their hard spectrum made it difficult to explain 
these events as shock breakouts, before it was realized that the radiation will be away from 
thermal equilibrium^. Plugging 10 48 erg and 150 keV, the total energy and typical frequency 



7 For example, iTan. Matzner fc McKed (|200ll ). which consider shock breakout as the source of energy 
of GRB 980425, add an interaction of the expanding gas with external mass shell, at a radius where the 
gas is optica ll y thin , in order to reprocess the gas energy into internal energy and radiate it non-thermally. 



Wang et al. I (|2007f ) suggested that the gamma-rays are produced by bulk comptonization of low energy 
thermal photon s during the breakout. A proce ss that does not work at the relevant temperatures due to 
Compton loses (jKatz. Budnik fc Waxmanll2010j ) 
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of GRB 980425, into equation [18] results in a breakout duration of ~ 10 s, comparable 
to the observed one. The breakout radius in that case is ~ 6 x 10 12 cm and the final 
Lorentz factor i^f| jf >0 « 3. The associated supernova of GRB 980425, SN 1 998bw, was very 



energ etic exhibiting M e j « 15M and a kinetic energy of ~ 0.5 x 10 53 erg ( llwamoto et al. 



19981 ). Interestingly, Setting i?* ~ 6 x 10 12 and M ej = 15M in equation [26] we find that 
an explosion energy of E f» 3 x 10 53 produces a breakout with the observed properties of 
GRB 980425. This energy is slightly larger than the one seen in the SN 1998bw, but it is 
close enough to suggest that the energy source of the SN explosion is also the one of the 
shock breakout. The mild discrepancy between the two energies may be a result of the 
actual stellar density profile (compared to the power-law with n = 3) or of a deviation of the 
explosion from sphericity. Using the energy and typical frequency of GRB 031203, 5 x 10 49 
erg and > 200 keV, equation [TS] predicts a duration < 35 s. The observed duration, ~ 30 s 
is consistent with this value and indicates that the actual temperature in not much higher 
than 200 keV. The breakout radius in that case is ~ 2 x 10 13 cm and the final Lorentz factor 
is 7/ 5 o ~ 5. The SN associated with GRB 31203, SN 20031w, is similar in ejected mass 
and slightly more energetic than SN 1998bw ( IMazzali et al. 1120061 ). Setting i?» ~2x 10 13 



and M e j = 15M in equation [2^] we find that an explosion energy of E w 10 54 produces a 
breakout with the observed properties of GRB 031203. Again this value is larger by about 
an order of magnitude compared to the one observe in the associated SN. 

The breakout radii that we find for GRBs 980425 and 031203 are larger than those of 
typical Wolf-Rayet stars, the probable progenitors of GRB-SNe. This implies that either 
the breakout is not from the edge of the progenitor star (e.g., from a wind) or that low- 
luminosity GRB progenitors are larger than commonly thought (which may help chocking 
the relativistic jet). Finally, a prediction of relativistic shock breakouts is a significant X-ray 
emission, that can be comparable in energy to that of the breakout 7-ray flare. The X-rays 
are emitted during the spherical phase, over time scale that is significantly longer than that 
of the breakout emission. GRB 031203 has shown a bright signal of dust scattered X-rays, 
which indicates on ~ 2 keV emissio n with energy comparable to that of th e gamma-rays, 



during the first 1000 s after the burst ( jVaughan et al. 1120041 ; iFeng &: Foxll2010f ). These values 



fit the predictions of a breakout model. For example, if the breakout is from a stellar surface 
at = 2 x 10 13 and M ej - = 15M Q , then equations I3"U][3"21 predict a release of ~ 10 50 erg at 
~ 2 keV within the first hour after the burst. 



8 Although this value of the Lorentz factor is below the lower-limit derived bv lLithwick fc Sari (2001) for 
GRB 980425, there is no contradiction between the two studies. The reason is that iLithwick fc Saril (|200lh 
derived their limits assuming that there is a hard spectral power-law that extends to very high energies, 
similarly to the one observed in long GRBs. However such power-law does not exist in the case or relativistic 
breakout and therefore their constraints are not applicable to the model that we discuss here. 
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To conclude, we find that the energy, temperature, and time scales of all four low- 
luminosity GRBs are explained very well by shock breakout emission and that all of them 
satisfy the relativistic breakout relation between E bo , T bo and tf^ (equation [TS]). These 
observables are largely independent of the pre-shocked density profile and are therefore ap- 
plicable also for a breakout from a shell of mass ejected by the progenitor prior to the 
explosion and possibly to a breakout from a wind. Breakout emission can also explain many 
other observed properties, such as the smooth light curve, the afterglow radio emission, the 
small fraction of the explosion energy carried by the prompt high energy emission, and the 
delayed bright X-ray emission seen in GRB 031203. We therefore find it to be very likely 
that indeed all low-luminosity GRBs are relativistic shock breakouts. 



5.5. GRB 101225A 



GRB 101225A is a peculiar explosion showing an ~hour long smooth burst of gamma- 



: t erglow , which is consistent with a blackbody 
(201 lj) a rgue for a cosmo logical origin and 



rays and X-rays followed by a pec uliar IR-UV a '. 
emission. Its origin is unclear. iThone et 

estimate its redshift to be ~ 0.3, while iLevan &: Tanvirl (120111 ) argue that its origin is local. 
Here we do not attempt to determine which, if any, of the two views is correct. Instead, 
given the similarity of the high energy emission of GRB 101225A to some low-luminosity 
GRBs, we ask whether shock breakout can be the sour ce of this emis s ion, a ssuming that the 
burst is cosmological and that the redshift estimate of IThone et al. I (120 111 ) is correct. 



We test if its observed properties satisfy equation [18j which with a total energy of ~ 10 51 
erg and gamma-ray spectrum that peaks at ~ 40 keV predicts a duration of ~ 10 4 s. This 



value is consistent with the constraints on the actual duration of th e flare (IThone et al. 
201 ll ). The resulting breakout radius is ~ 3 x 10 14 cm. Interestingly IThone et al. I (120111 ) 



present a model that attempts to explain the IR-UV afterglow by an expanding opaque shell 
with velocities at least as high as ~ 1/3 and initial radius of 2 x 10 14 cm. A breakout 
of a shock that accelerates such a shell pro duces a high energy signal that is similar to the 
observed one. We therefore conclude that if IThone et al. I ( 120111 ) redshift estimate is correct, 
then a relativistic shock breakout probably provides the simplest explanation for the high 
energy emission. 
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5.6. Long gamma-ray bursts 



In the colapsar model of long GRBs the relativistic jets pierce through the envelope 
and emerge in order to produce the observed gamma-ray emission. During their propa- 
gation the jets drive mildly relativistic bow shocks into the envelope and inflate a cocoon 
that colli mates the jets to an even na r rower angle th a n their initial launching opening an- 



gle (e.g., iMacFadven fc Woosley 



1999 



Matzner 



2003 



Morsony. Lazzati fc Begelmanl 12007 



Mizuta fc Aloyl 120091 ) . Recently, iBromberg et al. I (120111 ) derived an analytic description of 
this interaction while the jet is propagating deep in the envelope, before it encounters the 
sharp density gradient near the stellar edge. 

The term shock breakouts from GRBs is used in the literature to describe several dif- 
ferent phenomena. Here, similarly to the rest of the paper, we focus on the emission that 
is escaping from the breakout layer after the forward shock, which is driven by the jet head 
and the cocoon into the stellar envelope, breaks out. The structure of the shock that is 
driven into the stellar envelope is highly non spherical but at the time of the size of causally 
connected regions is small an a local spherical approximati on is appropriate. Thus , the the- 
ory discussed above, together with the mod el presented in 
results seen in numerical simulations (e.g., iMizuta fc Alo\ 
rough idea of the expected breakout signal. 

We consider a typical GRB jet with a half-opening angle 8q = 0.1 rad and a total 
luminosity of 10 50 erg/s. The jet is launched continuously into a 5M progenitor with a 
radius of 5R Q . At the point that the jet reaches the steep density drop (i.e. r ~ R*/2) 
i ts velocity approaches the speed of light and it is narrowly collimated to an angle of ~ 6q 



. 3romberg et al. I (120111 ) and the 



20091 ). can be used to obtain a 



( IBromberg et al. 1 1201 lh . The opening angle of the cocoon at this point is ~ 8 . From 
that point the jet head and the cocoon shock starts accelerating towards the stellar edge 
and to expand sideways. Numerical simulations show that the head is expanding sideways, 
becom ing comparable in siz e to the cocoon which does not spread significantly (e.g., figure 
17 in IMizuta fc Aloyl 120091 ). Once the fractional distance of the shock to the edge, z/R*, 
is comparable to the shock opening angle, in our case ~ 6q = 0.1, then the shock lose 
causal connection with t he sides and it can be ap proximated by the spherical solution of an 
accelerating blast wave ( IJohnson fc McKedll97ll ). Taking r y s f3 s (z = O.li?*) ~ 1 and using 
the relation 7 S /3 S oc z~ om with equation we obtain 7 ~ 10. This value is larger than 
7o jS so the breakout shell accelerates also after the transition to the spherical phase. We 
do not provide here the solution to such case, but the initial energy in the breakout shell 
is ~ lO 48 (-R*/5-R ) 2 erg and its energy increases during acceleration by at most an order 
of magnitude. Assuming that there is no significant wind surrounding the progenitor, this 
energy is released by a short, ~ ms, pulse of >MeV photons. 
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The main distinguishable feature of the breakout pulse is its harder spectrum and low 
fluence compared to the total burst emission. The pulse luminosity depends strongly on 
the progenitor size. If the progenitor is rather large, > 5i? , then the breakout can be 
detected by Swift and Fermi up to a redshift ~ 0.1. If it is much smaller, then the breakout 
is too faint for detection. Even in cases that the breakout is detectable, it may be difficult 
to separate it from the GRB itself. Nevertheless, it is worth looking for an initial hard spike 
in nearby GRBs. Interestingly, there is a su b-sample of BATSE GR Bs that show an initial 
short hard spikes followed by a softer burst ( jNorris fc Bonnelll 120061 ) . 



6. Summary 

We calculate the emission from mildly and ultra relativistic spherical shock breakouts. 
First, we determine the pre and post breakout hydrodynamical evolution. This requires 
knowing the shock width, which we obtain by deriving an improved analytic description to 
the struc ture of the transition layer in relativistic radiation mediated shocks, based on the 



results of iBudnik et al. I (120101 ) . After obtaining the hydrodynamical evolution we follow the 
temperature and opacity evolution to find the observed light curve. Our calculations are 
applicable to stellar explosions with shock velocity larger than 0.5c, in which the breakout 
shell ends its acceleration during the planar phase, i.e., its final Lorentz factor < 30 for 
typical stars. We consider only cases where the progenitor wind has no effect on the observed 
emission. 

A relativistic radiation mediated shock bri ngs the gas to a roughly constant rest frame 



temp erature, ~ 200 keV, and loads it with pairs (IKatz. Budnik &: Waxmanll2010uBudnik et al 



20101 ). Photons remain confined to the post-shock expanding gas as long as pairs keep it 
optically thick. A significant number of photons are released towards the observer only from 
optically thin regions, which has a gas rest frame temperature < 50 keV (low pair load) and 
unloaded pair gas opacity, r, smaller than unity. Our solution of the shock structure shows 
that the shock width is significantly smaller than r = 1, implying that the shock accelerates 
also shells with r < 1, which become transparent once their pairs annihilate. Among these 
shells the one that carries most of the energy is the most massive one, i.e., r 1, which we 
refer to as the breakout shell. In the regime that we consider the breakout shell becomes 
transparent during the planar phase (i.e., before its radius doubles) and after it achieves its 
final Lorentz factor. Light travel time and relativistic effects dictate that the emission of 
the breakout shell, once it become transparent, dominates during the entire planar phase. 
Slower shells release their photons only during the spherical phase. The processes described 
above produce the main following observables: 
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(i) A j-ray flare with a typical photon energy that ranges between ~ 50 keV and ~ 2 MeK 
TTie /Zare typically contains a small fraction of the explosion energy and it may be signifi- 
cantly shorter than the progenitor light crossing time. This flare is generated by the breakout 
shell when it becomes transparent. 

The Energy, temperature and duration of the flare depend only on two parameters, 7^0 and 
the breakout radius. Thus, in addition to providing the value of these two physical parame- 
ters, the three observables must satisfy the relativistic breakout relation (equation ITS]) . This 
relation can be used to test if an observed flare may have been produced by a quasi-spherical 
relativistic breakout. Note, that our calculation of the breakout flare energy, temperature 
and duration are independent of the exact profile of the pre-breakout density, as long as it 
is steeply decreasing. 

(ii) Due to light travel time effects, the breakout flare contains photons from faster and 
lighter shells, producing a high energy power-law spectrum, vF v oc z/ -0 74 . It also contains 
photons emitted by the breakout shell after it becomes transparent and cools adiabatically, 
producing a low-energy power-law spectrum, vF v oc v. Both power-laws are limited to about 
one order of magnitude in frequency. 

(iii) The flare ends with a sharp decay at a time that coincides with the transition to the 
spherical phase. If the breakout is ultra relativistic (i.e., 7^0 3> 1) then L oc t~j? s and T oc t~^ s , 
until the spherical phase emission dominates. 

(iv) After the flare ends, spherical evolution dictates a steady decay of the luminosity. At 
first, while relativistic shells dominate the emission, L oc t~ bs l . Later, at t$f^, once New- 
tonian shells become transparent L oc t^' 35 . Note that during the latter phase the total 
emitted energy increases with time. 

(v) The post flare temperature decays at first at a steady rate of about £~ fe °' 6 , generated 
during the spherical phase. A sharp drop in the temperature is observed when shells which 
were not loaded by pairs during the shock crossing (v s < 0.5) dominate the emission. The 
drop is typically from the X-ray range to the UV, and it is observed at about 10 t ^. 

(vi) The energy emitted before the steep temperature drop is often comparable to that emit- 
ted during the breakout flare. Thus, a signature of relativistic breakout in many scenarios is 
a bright and short / y-ray flare and a delayed X-ray emission with comparable energy. 

We apply our model to a range of observed and hypothesized explosions finding the 
following predictions: 

• Type la SNe and other white dwarf explosions: A range of white dwarf explo- 
sions, including observed la SNe and SNe 2005E & 2010X as well as hypothesized .la 
SNe and AIC, are all predicted to produce a similar breakout emission. A ~ 1 — 30 
ms pulse of 10 40 — 10 42 erg composed of ~ MeV photons. Such pulse from a Galactic 
explosion is detectable, implying a low, yet not negligible, chance to detect such flare 
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from type la SNe in the coming decade. The rate of the other types of explosions are 
not known. Moreover there is almost no correlation between the properties of the 7-ray 
flare and the luminosity of the optical emission. Implying that Galactic breakouts from 
white dwarf explosions may have already been observed, but not recognized as such. 
In this context the suggested sub-class of very short GRBs (ICline et al. 1 120051 ) may 
contain such events. 

Broad-line Ib/c SNe: Their progenitor is compact enough and the energy of some is 
high enough to form mildly relativistic breakouts. For example SN2002ap-like events 
are predicted to produce a detectable breakout 7-ray signal in case that there is no 
thick wind surrounding the progenitor. 



Extremely energetic SNe: SNe such as 2007bi ( IGal-Yam et al. 1120091 ) and explo- 
sions of the type reported by iQuimby et al. I ( 120091 ). are energetic enough to generate 
mildly relativistic breakouts in very massive stars. These breakouts are likely to release 
10 44 — 10 46 erg in soft 7-rays over seconds to minutes. Such flashes are detectable in 
the nearby Universe, however these are rare events, and the chance to have one in a 
detectable distance is low. 



Low luminosity GRBs: The origin of low-luminosity GRBs is unknown. iBromberg. Nakar fc Piran 
( 120111 ) show that low-luminosity GRBs are very unlikely to be produced by relativistic 
jets that pierce through their host star, which is a necessary ingredient of the collap- 
sar model for long GRBs. In that case what can be the source of the high energy 
emission? The smooth light curve and the mildly r elativ istic ejecta that generate the 
radio emission of GRB 980425 lead lKulkarni et al. I (119981 ) and lTan. Matzner fc McKee 
( 120011 ) to suggest that relativistic shock breakouts are rela ted to low luminosity G RBs. 
The detection of a thermal component in GRB 060218 (jCampana et al. 1 120061 ) and 
of additional low luminosity GRBs with similar proper t ies to GRB 980425 , supported 
this suggestion (IWaxman. Meszaros &: Campanal 120071 : IWang et al. 1120071 ). especially 
after i t was realized that the shock em ission strongly deviates from thermal equilib- 
rium ( IKatz. Budnik fc Waxmanl 120101 ). However, there was no quantitative model 
of the breakout emission that includes deviation from thermal equilibrium and gas- 
radiation relativistic dynamics, which could test this suggestion. Moreover, it was 
unclear whether shock breakouts can generate the energies seen in these bursts (espe- 
cially GRBs 980425 and 031203 due to their shorter duration). 

We find that relativistic shock breakouts can explain very well the energy, temperature 
and time scales of the prompt emission of all observed low-luminosity GRBs. Moreover, 
all four known bursts satisfy the relativistic breakout relation between these three 
observables (equation [TS]) . We also show that relativistic breakouts provide a natural 



-28 - 



explanation to the low gamma-ray luminosity compared to the total explosion energy 
and to the hard to soft spectral evolution (including the late energetic X-ray emission 
that inferred from dust echoes of GRB 031203). We therefore find shock breakouts to be 
the most promising sources of all low-luminosity GRBs. The connection between low- 
luminosity and long GRBs (both share a relation to a similar type of SNe), suggests 
that the energy source for the shock breakout are choked relativistic jets that are 
launched by a GRB engine but fail to punch through the star. 

In this context we examine the high energy emission of GRB 10 1225 A, whose origin 
(local or cosmo logical) is undetermined yet, but its prompt e mission resemble low - 



luminosity GRBs. We find that if the redshift determination of iThone et al. I ( 120111 ). 
z ~ 0.3, is correct then this burst also successfully passes the test set by equation 
[T8| suggesting that the source of its high energy emission is a mildly relativistic shock 
breakout at a relatively large radius of ~ 3 x 10 14 cm. 

• Long GRBs: According to the popular view of long GRBs, relativistic jets drill their 
way through the progenitor envelope driving a forward shock into it. We estimate that 
a typical value of the initial breakout Lorentz factor of this shock is 70 ~ 10. Our 
theory cannot determine its final Lorentz factor, since it ends its acceleration after 
the evolution enters the spherical regime. Nevertheless, this value of 70 indicates that 
shock breakout from long GRBs release ~ lO 48 (-R*/5-R ) 2 erg in a short pulse of >MeV 
photon. This pulse is too dim to be observed in the typical redshift of long GRBs, but 
it may be observable in nearby z ~ 0.1 GRBs. 

We conclude that relativistic shock breakouts are a generic process for the production 
of gamma-ray flares, which opens a new window for the study and detection of a variety of 
stellar explosions. 

We thank Omer Bromberg, Amir Levinson and Tsvi Piran for helpful discussions. E.N. 
was partially supported by the Israel Science Foundation (grant No. 174/08) and by an IRG 
grant. R.S. was partially supported by ERC and IRG grants, and a Packard Fellowship. 



APPENDIX - The transition layer structure of relativistic radiation mediated 

shocks 



Here we derive an analytic descripti on of the tran s ition layer of relativistic radiation 
mediated shocks, b ased on the re s ults o f iBudnik et al. I ( I2010l ). We use mostly similar no- 
tations to those of IBudnik et al. I (120101 ) and carry out the calculation in the shock frame, 
where the flow is assumed to be in a steady state. We approximate the gas-radiation in 
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the transition region as two counter flows. One is a flow of photons with energy ~ m e c 2 
that are generated in the immediate downstream and flow towards the upstream at a drift 
velocity ~ c. The other is a gas flow that is moving towards the downstream. This is a flow 
of the fluid that is shocked, which is moving with a Lorentz factor T u in the far upstream 
(this is the shock Lorentz factor) and is slowed down in the transition region to /3 = 1/3 
(i.e., T 1). The downstream moving fluid is composed of protons with a proper density 
n and pairs and photons with a proper density n± tT The photons in this flow were part of 
the photon flow that moves towards the upstream, before they were scattered once by the 
fluid leptons and joined the fluid motion towards the downstream^. We are interested in 
the shape of the transition, i.e., T and the pair fraction as a function of optical depth and 
physical length. Our final goal is to find the shock width in the upstream frame in units of 
the pre-shock (i.e., with no pairs) upstream Thomson optical depth. 

The shock is generated by the interaction of the two counter-streaming flows. Every 
collision between upstream moving photon and a downstream moving fluid lepton [photon] 
transfers energy and momentum to the photon [created lepton] and makes it joining the 
downstream moving fluid. The temperature of the pairs and fluid photons is ~ Tm e c 2 , as 
seen in the fluid rest frame. In a steady state the flux of photons that are moving towards the 
upstream, n 7 -+ us , is equal to the number flux of pairs and photons that are moving with the 
fluid (assuming that the later dominate over the proton number flux), i.e., n^^ us « r/3n± )T 
Now, since upstream moving photon that is scattered is joining the downstream moving 
fluid, d(Tn± )1 ) = —n 7 ^ us dr s , where t s is the optical depth for photons moving toward the 
upstream, defined to be increasing in the upstream directiorj^l. Defining x± ~ = n± t7 /n and 
assuming T 2 >> 1 (i.e., (3 ps 1), we obtain: 

^T 2 = -*±,r (A-l) 
dr s 

Number and momentum flux conservation imply: 

n u T u (5 u = nTP (A-2) 

r^A^m^c 2 = T 2 f3nc 2 (m p + Ax±^Tm e ) (A-3) 

In the momentum equation we assume that the momentum flux of the photon field that 
streams towards the upstream is negligible, which is true as long as Y 2 1. We also assume 



9 Protons and pairs are assumed to be coupled by collective plasma processes on scales much shorter than 
any scale on interest in the system and thereby can be treated as a single fluid. The photons in this flow are 
having the same downstream drift motion, but are not directly coupled to the protons and leptons. 

(|2010h 



10 



Our definition of r, is similar to — r in the notation of Budnik et al. 
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CO. 




Fig. 2. — Comparison of the shock structure for four values of T u (6, 10, 20 & 30), obtained 
by integr ating over equa t ion IA-61 fusing r* = t S (ji r J^ AT ^ ), [thin solid line], to the numerical 
results of iBudnik et al. I ( 120101 ) (their figure 6), [thick dashed line] 



x± rj ^> 1. These equations imply: 



x±. 



m p T u -Y 
4m e T 2 



Therefore: 



dT dT dx±, 
dr s dx± n dr s 



(A-4) 



(A-5) 



2T U - r 

If we define r* as the Thompson optical depth of a photon moving toward the upstream then 
where o~t is the Thomson cross-section and o~kn(^ 2 ) is the Klein- Nishina 



(7 T 



'* Is °kn{T 2 )' 

cross section of electron at rest and photon with energy T 2 m e c 2 . Therefore 



2r„ - r 



Ta KN (T 2 ) T u -T 



dT 



(A-6) 



Integrating this equation starting at the end of the transition layer, i.e., r(r s = 0) = 1 results 
in the structure of the t ransition layer. Comparison of this structure to the numerical results 
of IBudnik et al. I ( 120101 ) is presented in figure [2j 



In order to find the width of the shock in the upstream frame we use the relation 

fir* = Y u n u x ±1 a T dz' (A-7) 



fc* 0.5 
U 
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n g t z' (1 +ln[2r 2 ])/r 2 
u T v 1 u J/ u 



Fig. 3. — The structure of the shock transition region of relativistic radiation mediated 
shocks, presented in an optical depth scale (n u aT z ' (1+2 ln[r^])/r^), which makes it universal, 
i.e., the value of T/T u is independent of T u (for values larger than 1/F U ). 



where z' is the length coordinate (positive towards the upstream) in the shock frame. There- 
fore 



4m P 



2r u -r 



-rdr 



(A-8) 



m p (T u - T) 2 r u a KN (T 2 ) 

Integrating dT from the shock towards the upstream, the value of T/T u (for values larger than 
l/r u ) has a universal profile (independent of T u ) as a function of n u OTz'{\ + 21n[r^])/r 2 . 
This profile is presented in figure [31 This profile implies that the shock width in unit of pair 
unloaded Thomson optical depth in the shock frame is roughlj0, oc T 2 . The value of the 
proportionality coefficient is somewhat arbitrary and it depends on the value of T/T u that 
defines the "shock width". We chose r/r u = 0.9 obtaining n u aTz' ~ O.Oir 2 (the value of 
the logarithmic factor is 5 — 10 for Y u = 2 — 10). Choosing a different value of r/r„, to define 



11 Note that here we find the shock width in pre-shock, witho ut pairs, Thomson opt ical depth, which can be 
directly related to physical width for a given upstream density. iBudnik et al. I (|2010l ) find the shock width in 



units of total Thomson optical depth, including pairs, of a photon that crosses the transition layer towards 
the upstream, t* . As it turns out, in both cases the width oc T 2 but the proportionality coefficient is different 
(IBudnik et al. Il201ol find At* « T 2 ). 
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the shock width, changes the coefficient, which in turn has a very weak effect on equation [HJ 
Now, since the length scale in the upstream frame, z, is shorter by a factor T u than in the 
shock frame (i.e., z' = T u z), and defining t u = n u <7T z , we approximate the shock width in 
units of Thomson optical depth of the pre-shock upstream (i.e., with no pairs) as: 

Ar u ~ 0.0ir u (A-9) 

Note that our analytic description is applicable to shocks were ^> 1 and are therefore not 
directly applicable to mildly relativistic shocks. We do expect however that pair creation 
in these shocks, which increases the optical depth per proton in the immediate downstream 
by a factor of the order of m p /m e , will result in Ar u <C 1 and that equation IA-91 can be 
extrapolated to the mildly relativistic regime. 
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